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Abstract 

We analyze temporal and spatial meson correlators in quenched lattice QCD at T > 0. 
Below T c we observe little change in the meson properties as compared with T = 0. Above 
T c we observe new features: chiral symmetry restoration and signals of plasma formation, 
but also indication of persisting "mesonic" (metastable) states and different temporal and 
spatial "masses" in the mesonic channels. This suggests a complex picture of QGP in the 
region 1 — 1.5 T c . 
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1 Introduction 



With increasing temperature, hadronic correlators are expected to change their nature dras- 
tically (see, e.g., [0, [§])• At the critical temperature, the deconfinement of color degrees of 
freedom and the restoration of the chiral symmetry are expected to occur simultaneously. Two 
"extreme" pictures are frequently used to describe the low and the high T regimes, respectively: 
the weakly interacting meson gas, where we expect the mesons to become effective resonance 
modes with a small mass shift and width due to the interaction; and the perturbative quark 
gluon plasma (QGP), where the mesons should eventually disappear and (at very high T) 
perturbative effects should dominate. 

Near to the critical temperature, however, the actual physical situation is more involved. 
The interaction with a hot meson gas and with baryonic matter have been studied in various 
phenomenological models which predict appreciable changes in the vector meson properties (see, 
e.g., H). In a NJL model [f| the scalar and the pseudo-scalar modes are found to correlate 
strongly, and to subsist even above the transition as so-called "soft modes" , corresponding to 
narrow peaks in the spectral function and realized as the fluctuation of the order parameter of 
the chiral symmetry restoration transition. On the other hand, at the short distance scale, the 
fundamental excitation should be quarks and gluons. Lattice QCD results on the quark number 
susceptibility support this view [||]. These pictures may not be contradicting each other: DeTar 
conjectured the existence of excitations in the QGP phase corresponding to different distance 
scales, and pointed out to the possibility of "confinement" still ruling the large distance scales 
||. While with increasing temperature the physics should appear increasingly dominated by 
quark and gluon degrees of freedom, in accordance to the perturbative high temperature picture, 
the intermediate temperature range above T c seems to be much more complex and dominated 
by strongly interacting quarks (see also @), which, as we shall see below, even tend to stay 
strongly spatially correlated and thus agree with a picture of effective, low energy modes in 
the mesonic channels. Since these matters are related to questions about the evolution of the 
early universe, on one hand, and to the interesting results from heavy ion collision experiments 
||], where QGP conditions are being realized Q, on the other hand, it is important to have 
quantitative estimates in addition to qualitative understanding and we need model-independent 
studies of the hadronic correlators at finite temperature. 

Lattice simulations are the most powerful instrument at present to investigate such prob- 
lems in the fundamental theoretical framework of quantum chromodynamics (QCD). Exten- 
sive studies have been dedicated to the thermodynamics of the finite temperature transition 
(see, e.g., pO[|,pT[|). Concerning the hadronic sectors, numerical analysis of "screening" (spa- 
tial) propagators indicate correlated (bound?) quarks while the mesonic "screening masses" 
increase toward the 2-free quark threshold (2ttT, induced by the anti-periodic boundary con- 
ditions in the temporal direction) |fl2j l. Since the spatial directions can (and must) be made 
large, these propagators are unproblematic in principle and can be studied as well as for T = 0: 
at any T > the propagation in the space directions represents in fact a T = problem (with 
asymmetric finite size effects). The interpretation of the results from "screening" correlators 
in terms of modes of the temporal dynamics is far from straightforward, however: since in the 
Euclidean formulation the 0(4) symmetry is broken at T > (see, e.g., physics appears 

different depending on whether we probe the space ("<r": x) or time ("r": t) direction. The 
static quark potential associated to propagation in a spatial direction, for example, is a very 
anisotropic quantity which above T c still grows linearly in two of the spatial directions (con- 
fines), in contrast to the potential associated to the propagation in the t— direction which is 
isotropic and not confining. Therefore we need to investigate hadronic correlators with full 
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"space-time" structure, in particular the propagation in the Euclidean time. 

Ideally, what we should like to do is to reconstruct the spectral function in the given channel. 
Then we could directly compare with the results from heavy-ions experiments, see e. g., ||. The 
spectral function at finite temperature can be extracted from the correlators in the (Euclidean) 
temporal direction whose extent l T is related to the temperature as T = l/l T |l4|, [D|. These data 
(after Fourier transforming the t-correlators) are given at discrete Matsubara frequencies on the 
imaginary energy axis and are affected by errors. The extraction of the spectral function implies 
(logically) an interpolation and an analytic continuation to the real energy axis. For a numerical 
analysis, which produces a finite amount of information, this is an improperly posed problem. 
Its solution is dependent on imposing supplementary conditions ("a priori information") to 
regularize the algorithms and to prevent the amplification of the errors. These conditions can 
be either based on general, statistical arguments (e.g. variance limitation, Bayesian analysis, 
maximal entropy method) or on particular, physical expectations (e.g., using an ansatz for the 
spectral function which leads to an explicit analytic form for the correlator, to be fitted against 
the data). We should, however, be aware of the fact that all regularization introduces a bias 
and therefore this problem is fundamentally intricate. 

The main difficulty in the numerical calculation at T > originates in the short temporal 
extent l T = 1/T. Beyond the general necessity of producing enough and precise data the T > 
problem is doubly complicated as compared to the T = one: on one hand the structure 
we may expect is more complex than just a pole, on the other hand the time extent of the 
propagation cannot be made large to select the low energy contributions. We shall now briefly 
discuss this questions and thereby also introduce our procedure to deal with them. 

1) Lattice problems: Large T can be achieved using small N T = l T /a (a: lattice spac- 
ing), however this leads to systematic errors [16]. Moreover, having the i-propagators at only 
a few (N T /2) points makes it difficult to characterize the unknown structure in the corre- 
sponding channels: practically any ansatz can be fitted through 2-3 points! To obtain a fine 
t-discretization and thus detailed i-correlators, while avoiding prohibitively large lattices (we 
need large spatial size in order to avoid finite size effects, typically l a ~ 3l T ), we proposed 
[17, 18] to use different lattice spacings in space and in time, a a /a T = £ > 1 |19|. The renor- 
malization analysis of such lattices, however, is more involved, because of the supplementary 
parameter £ and this introduces also some uncertainties. 

2) Physical problems: The low energy structure of the mesonic channels cannot be observed 
directly, due to the inherently coarse resolution 2ttT = 2ir/l T of the imaginary energy axis. 
Refining the discretization of the time axis improves the fitting and analytic continuation 
problem, but although we are following the question of the spectral analysis we do not have yet 
reliable results at T > for this challenging question. Our problem setting here is therefore 
more limited: we shall try to recognize mesonic states and ask about their character and 
properties at various temperatures. 

In this paper, we investigate the full four-dimensional structure of the meson correlators 
on an anisotropic lattice in the quenched approximation. Thus the phase transition is the 
deconfining one, and the hadronic correlators are constructed with quark propagators on the 
background gauge field. 

Our strategy here is the following: we first select the mesonic ground states of the T = 
problem (where the time direction can be made sufficiently large - about 3.2 fm in our case, 
which, at the quark masses we work with, means about 8 pion correlation lengths) and charac- 
terize their internal structure by measuring the (Coulomb gauge) wave functions. Then we ask 
whether states characterized by similar internal structure can be retrieved at higher temper- 
ature, try to reconstruct them with help of correspondingly smeared sources and investigate 
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how they are affected by the temperature. If the changes in the correlators are small, which is 
consistent with mesons interacting weakly with other hadron-like modes in the thermal bath, 
this procedure allows us to define "effective modes". Large changes will signal the breakdown 
of this weakly interacting gas picture and there we must try to compare our observations with 
other pictures, in particular the perturbative QGP. 

The meson correlators in the temporal direction play a central role in this study, which 
is therefore meant to supplement other approaches, including studies of screening propagators 
[12], [13]. To understand the effect of fixing a mesonic source we employ three kind of meson 



operator smearing. The propagators and the wave functions are also compared with those of 
mesons composed of free quark propagators ("free" mesons) []. Finally we attempt a chiral limit; 
note, however, that even with anisotropic lattices the short physical extent in the temporal 
direction makes the quantitative estimate of the (temporal) masses (if they exist) difficult. 
Our program should not be understood as an alternative for a study of the spectral function 
at T > 0, but as an attempt to answer some special questions about the phenomena in QGP. 
In that sense our results only offer partial views. 

To prevent a certain confusion we stress here that we do not look for the eigenstates of the 
Hamiltonian (transfer matrix), which show up as asymptotic states for t — > oo at T = 0. At 
non-zero temperature the physical processes are essentially dependent on the mixtures induced 
by the thermal bath. For the physical picture and for building models the question is whether 
these phenomena can be described in terms of some effective excitations (quasi-particles ]pO|l ), 
which "replace" thus the fundamental particle modes, or completely new states dominate the 
physics above some T. 

Our analysis proceeds in three steps: 

1. Analysis of temporal propagators. Here we try to see what kind of excitations propagate 
in the mesonic channels at T > based on the t— dependence of these propagators 
("effective mass" m{t)). 

2. Analysis of the "Coulomb gauge wave functions". Here we study the behavior of the 
temporal correlators with the distance between q and q at the sink, which provides us 
with information of the spatial correlation between the quarks at given t. 

3. Analysis of the temperature dependence of the temporal and spatial masses of the putative 
states which are compatible with the behavior observed at the previous steps. 

Note that due to the quenched approximation the dynamics is incomplete. In a strong 
sense the Hamiltonian does not posses true mesonic states and only provides the gluonic inter- 
actions responsible for the forces binding the quarks. This is not a problem specific to nonzero 
temperature but is the same already at T = 0. The success of the spectrum calculations at 
T = indicates, however, that one should not consider quenched QCD as a theory by itself but 
as an approximation to the full theory (which posses genuine asymptotic mesonic states) and 
the exclusion of q-q pair creation as a reasonably small error at least concerning some of the 
characteristics of the hadrons. In particular, we observe strong indication of chiral symmetry 
restoration above the transition temperature. 

The paper is organized as follows. In the next section we describe our analysis strategy 
in some more detail. Section || describes the preparation of the lattice: the introduction of 
anisotropic lattice actions and the simulation parameters (the "calibration", i.e. the tuning 
of the anisotropy parameters, is described in detail in the Appendix). The subsequent two 

1 These can be seen as quark-antiquark correlation functions in the corresponding meson channels in lowest 
order perturbation theory. 
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sections present the results of the simulation: In Section 4 we observe the correlators at zero 
temperature and discuss the source smearing and the variational analysis. The results at finite 
temperature are presented in Section [|. The last section is reserved for discussion and outlook. 



2 Analysis strategy 

2.1 Comments on the physical problems 

We here should like to illustrate the problems raised by the finite temperature and the question 
of the source in the frame of our approach. The reader who is familiar with these problems 
may skip this section. 

Let us consider that we use some meson operator <£, then the propagator at j3 = 1/T < oo 
in Euclidean time t > is: 
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where (f) n are eigenstates of the Hamiltonian representing, say, (multi-)meson states and other 
hadron-like modes and we have: 



(0n|$|0fc) = C n k, C kn = C* nk = C nk (5) 

Here we expressed everything in units of a r (a~ 1 ), hence (3 = N T ; T = ex.p(—H) is the transfer 
matrix. For T = only the vacuum n = survives in the sum over n in (|3|). Assume each <I> 
selects not only the mesonic ground state, say |</>i), but also some other, excited states, then 

(0 fe |$|0o) ~ c i<5fc,i + c 2<5fc,2 + ■ • • (6) 

and the zero temperature propagator is (we put Eq = for simplicity): 

G^(t)= i cg 1 e-^+cg 2 e-** + ... (7) 

Hence the lightest state contribution will dominate at large t. Tuning a "perfect" source at 
T = we ideally achieve CQ k = for k ^ 1 and thus see only this contribution at all t. Suppose 
that we have been able to construct in this way a "perfect" operator <I>i ~ a + a) , with a(a^) 
the annihilation (creation) operator for a meson in the ground state. Then at T = G reduces 
to the first term, as desired 

G^-^)(i)^ ie -^ (8) 

(note that correlators with different operators at the source and the sink also project only on 
the ground state if either the source or the sink is "perfect"). With increasing T, however, 
further states beyond the vacuum survive in the sum over n and acting on each of them $i 
"adds/subtracts" a meson to whatever is there, correspondingly selecting from the inner sum 
the states k onto which this new state projects, (^nl^il^fc) 7^ 0, in a sloppy notation k £ {n±i}. 
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Instead of (|8|) the correlator is now a sum of contributions and we ask whether this mixture 
can be described by an effective mode of energy e[^ such that we can write, similarly to 
the T = expression: 



n ke{n ±1 } 







-PE„ 



cosh 



--t\E 



(9) 



To fix the ideas let us consider an oscillator with frequency u and a small anharmonic 
perturbation: 

H = H + eHj (10) 

(this may be considered a caricature of a weakly interacting meson gas, say). To first order in 
e we can use the unperturbed basis to calculate the propagator G(t). Let <3?i be the ground 
state operator, 



$i = (a + a^/y/fa , (0 fc |*J|0 n ) = -= (Vn+18 kn+1 + y/n6 k ,n-i) , 

\J Luj v ' 

then we have: 
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For the unperturbed oscillator (e = 0): 
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= E n — E n ^\ — uj. 



(13) 



Then the t— dependence factorizes in (12) and we have a trivial effect of the temperature: 
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If we turn on the interaction the levels are no longer equidistant (the effect of adding one more 
meson depends on the total number of mesons present in the state) and the t— dependence is 
non-trivially affected by temperature. We write: 



Ai = u, A n > 2 =u — e\ n , 
then to first order in e (weakly interacting gas) 



P i \ - 

t] UJ 
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Notice that the above effects show up although we use the "perfect" source $i: they represent 
the genuine temperature effects for an interacting system. From (|il^),(|l7|) we see that as long 
as the interaction between the modes ("the mesons") is weak we expect small changes which 
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may be simulated by a shift (and possibly a widening) of the peak in the spectral function, 
defining in this way an effective mode (^). Large changes, on the other hand, will signal the 
installation of a new regime. Then we must try to obtain additional information by other 
tests. Essentially, this is our program. Of course in real life we shall not be able to obtain a 
"perfect" source in the above sense. The various uncertainties inherent in our procedure will 
be repeatedly discussed in the course of the paper. 

If we use a "perfect" source but a different sink (with non-zero projection on the source) we 
reach similar expressions. To the next order in e, however, at T > the temperature correction 
to the mass will depend on the sink operator. Generally therefore at T > we expect to find 
sink dependence of propagators even for "perfect" source. This dependence can be seen as an 
indication for the importance of temperature effects. 

2.2 Mesonic correlators 

A first attempt to optimize the mesonic operators, in the spirit described in the previous section, 
is to introduce a smearing function u>(y), such that the zero-momentum mesonic operator reads: 

$tl\t)=^^uj(y)q(z,t) lM q(z + y,t). (18) 

z y 

giving rise to smeared correlators (we shall omit the index "/9" in the following): 



G^'\t) = <T> W(0)j> (19) 

= Eii,i2,3/"i,2/S u ' GHM2/2) ( Tr S{z 2 , 0; z x , t)^ya b S\z 2 + y 2 ,0; z[ + */i,*)757m 

Here (.) means summation over Yang-Mills configurations. S is the quark propagator and 
7m = {75, 71, 1, 7175} for M={Ps,V, S,A} (pseudo-scalar, vector, scalar and axial-vector, re- 
spectively). In the scalar sector, only the connected part of the correlator is evaluated. The 
Coulomb gauge is used to produce the quark propagators S - this is of course irrelevant for 
the 2/1 = 2/2 = expectation values. Note that generally we keep different smearing functions 
at the source and the sink. This will allow, given a certain basis of operators {$ a }, to perform 
a variational analysis in order to attempt a further optimization of the mesonic operator - for 
details see section 4. 

As smearing functions we shall use two different kinds, 

i) A "point" source (sink): 

u s {y) oc S(x- y). (20) 
This will be mainly used to study the x dependence of the correlator 

G$ (x, t) = ^ (t) = £ (q(z + x, t) lM q(z, t) ^ (0)) . (21) 

z 

at fixed t. Gyp(x,t) can be interpreted as "Coulomb gauge wave-function", it indicates 
the spatial correlation between quark and anti-quark. 

ii) The convolution: 

Uabiv) = Y.v^aiv)^^ + y) (22) 
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which is equivalent to using smeared quark and antiquark fields with smearing functions 
u>b and uj a respectively. We use here three kinds of quark smearing functions: 



io p (y) oc 5{y) ("point") (23) 
oj e (y) oc exp(-o|y| p ) ("ezp.") (24) 
u w (y) oc 1 ("wail") (25) 



In tuning the exponential source exp in (24) we shall use the parameters a,p from the 



observed dependence on x at large t of the temporal Ps wave function with point-point 



source at T ~ 0, G^\x,t) - see ElL The mesonic operator "exponentially" smeared 
both at the quark and the antiquark corresponds to a mesonic source in the relative q-q 
distance given by the convolution ([22]). Therefore the exp-exp smearing with a,p from the 
wave function implies a meson source typically wider than the measured wave function 
GSf ) (x,t),t»l. 

The t dependence of the temporal propagators G^~* u ^ (t) depends on the spectral functions. 
On a periodic lattice the contribution of a pole in the mesonic spectral function to the t- 
propagator is oc cosh[M(i — N T /2)] (this M is therefore called "pole- mass" ).0 A broad structure 
or the admixture of excited states leads to a superposition of such terms. Fitting a given t- 
propagator by cosh(m(t)(t — N T /2)) at pairs of points t,t + 1 

t) 

(28) 



G(t) 


cosh 


™ { 3 f (t) (N T /2-t)_ 




G(t + 1) 


cosh 


'rr$ f (t) (N T /2-t-l) 



defines an "effective mass" m e ff(t), which is a constant if the spectral function has only one, 
narrow peak. The effective mass is a rather sensitive observable which shows effects of the 
source dependence, widening of the spectral function or existence of excited states, without, 
however, allowing to differentiate among them. To the extent that the effective mass reaches 
a plateau and permits to define a "temporal mass" at T > 0, the latter connects directly 
to the (pole) mass of the mesons below T c , while above T c it will presumably help analyze 
the dominant low energy structure in the frame of our strategy. By contrast, using spatial 
propagators we shall extract the "screening mass" uS u \ 

Errors are estimated by the single elimination Jackknife method, unless otherwise notified. 
For details see section ||. For various comparisons we shall also use in (]i~9l ) , instead of the quark 
propagators S measured in each MC configurations, free quark propagators, defining in this 
way "free" mesons. 



3 Lattice setup 

In this section we describe the preparation of the lattice on which mesonic correlators at zero 
and finite temperature are calculated. 

2 More precisely, the relation between the slope parameter in cosh, say M and the position of the pole, M is 

M = £^2(cosh(M/£) - l) (26) 

(and correspondingly 

M = ln(M/£ + 1) (27) 

for fermionic propagators) ]l7|] . In the following we shall neglect these corrections, since they remain below the 
other uncertainties of our data. 
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3.1 Anisotropic lattice 

We use anisotropic lattices on which the spatial and the temporal lattice spacings are different: 
a a ^ a T [19]. The simplest generalization of commonly used Wilson actions for gauge and 
quark fields is obtained as follows. 
For the gauge field action, 

S G = — E f 1 - ^ ReTr u a( x )) + /^7G E f 1 - l ReTr U i*( x )) ( 29 ) 

7G ^<i<3 V 3 J x,i<3 V 3 J 

where 

U^ix) = U^{x)U v {x + p,)Ul(x + O)Ul(x) (30) 

The bare anisotropy parameter jg controls the ratio of spatial and temporal lattice spacings. 
For the fermion, 

SF = -^- y £q(x)K(x,y)q{y), (31) 
iK ° x,y 



»=i 

- k g lF [(1 - 74 ) U A {x) 5 x+ly + (1 + 74 ) U A {x - 4) S x _ ly ] , (32) 

where K a is the spatial hopping parameter and 7_p is the bare anisotropy for fermions]^] For 
later convenience, we define k as 

- = --2(7j.-l)=2(mo+4), (33) 

where mo is the bare quark mass parameter in units of the spatial lattice spacing. At a later 
stage, we shall carry out the chiral extrapolation in this 1/k. 

The actual anisotropy £ is defined using certain correlators, F, containing gauge and fermion 



fields. In general, 7g and are different from £ because of the interaction 21]. To obtain 
the desired value of £, one needs to tune the values of 7g and 7^ by requiring the isotropy of 
correlators in physical units: 

F a (z) = F T (t = £z), (34) 

where t and z are understood in the corresponding lattice units, a T and a a respectively. This 
renormalization procedure is called "calibration". Since we compare the temporal and the 
screening masses, it is important to obtain £ as precisely as possible, and verify that these two 
kinds of mass coincide at T = 0. 

In the case of dynamical quarks, a precise calibration requires a large effort, since we 
generally have four bare parameters ((3, 7c K a and 7^) and only three physical parameters £, 
a a and m q which can be varied, therefore the condition of physical isotropy implies a non-trivial 
constraint among the bare parameters pq ]. In a quenched simulation, however, the situation 
is much simpler. Here one generates gauge field configurations using certain values of (5 and 
7g and reads £ from gluonic quantities. Then the fermionic parameters are determined such 
that they give the desired quark (or hadron) mass and the same £ is obtained from hadronic 
correlators. 



3 Note that the Wilson term corresponding to this ansatz has not a Lorentz invariant naive continuum limit. 
Since this is an irrelevant operator this feature should not affect our results, in fact this ansatz is more efficient 
in damping the additional fermionic modes. This may be different for quantities where the Wilson term acts 
marginally. 
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N, 




0.081 4.05(2) 
0.084 3.89(2) 
0.086 3.78(2) 



20 0.1601 0.0389 0.1777( 8) 0.1962(10) 
20 0.1633 0.0276 0.1493( 9) 0.1747(12) 
30 0.1648 0.0223 0.1341(10) 0.1644(13) 



Table 1: Quark parameters, "jp is determined by the calibration using first N con f configurations. 
The error of k and m q is not estimated {p/p is fixed in successive calculation). The tabulated 
meson masses are the values obtained with smeared correlators described in the next section. 
All masses are given in units of a~ . 

3.2 Lattice parameters 

We use lattices of sizes 12 3 x N T , where N T = 72 (T = 0), 20 (below T c ), 16 and 12 (above 
T c ), with couplings (3 = 5.68, = 4.0, in the quenched approximation. ^ Configurations 
are generated with the pseudo-heat bath algorithm with 20000 thermalization sweeps, the 
configurations being separated by 2000 sweeps. In most cases, 60 configurations are used, 
except for the calibration. The gauge field is fixed to the Coulomb gauge. The calibration 
is described in detail in the Appendix. From the calibration of the gauge configurations we 
obtain a renormalized anisotropy £ = 5.3(1). The lattice cutoffs determined from the heavy 
quark potential are a' 1 = 0.85(3) GeV (a CT ~ 0.24 fm) and a" 1 = 4.5(2) GeV (a r ~ 0.045 fm). 

Table p] summarizes the quark parameters and gives the meson masses as determined in 
the next section. As a guide for the mass range we are concerned with, the quark mass 
is estimated by a naive relation m q = (k _1 — n~ )/2 using the critical hopping parameter 
k c = 0.17144(11). This simulation deals therefore with quarks in the strange quark mass region. 
The boundary conditions for quark fields are set to periodic and anti-periodic in spatial and 
temporal directions respectively, except for the calibration (at T = 0) where anti-periodic b.c. 
in all four directions are used. 

In the following we always display dimensionless quantities, that is, they are understood as 
given in the corresponding lattice spacings or their inverses. Since we have two lattice spacings, 
a a and a T , when we shall compare different quantities we shall use the translation a a = £a r , 
i.e., the numbers featured can be understood as given in units of a T (a^ 1 ). 

3.3 Finite temperature 

In |23|] we have measured the Polyakov loop susceptibility as a function of 7g for several values 
of N T . At N T = 18, the peak was found between 7 = 3.9 and 4.0, which means on our 7 = 4 
lattice that N t = 18 is very close to and just above T c . The estimated temperature for our 
values of N T are therefore T ~ 0, 0.93 T c , 1.15 T c and 1.5 T c for N T = 72, 20, 16 and 12 
respectively. 

We found that the configurations above T c stay in a single Polyakov loop sector during 
the whole updating history, and that the behavior of meson correlators strongly depends on 
the sector. The hadronic correlators feel the deconfining transition if they are taken in the 
real sector, but they appear not to "notice" it if they are taken on configurations in one of 
the complex sectors [^2j . With dynamical quarks, the Z% center symmetry is explicitly broken 
and the Polyakov loop prefers to stay on the real axis. Since we regard the quenched lattice 

4 The lattice described in this paper corresponds to Set-B data in our earlier reports El . 
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as an approximation to the dynamical lattice, we restrict our simulation to the case with real 
Polyakov loop sector. 



4 Zero temperature analysis. 

In the context of our strategy the analysis of the zero temperature correlators serves as foun- 
dation for the analysis at T > 0. This is also a good opportunity to describe in detail our 
procedure. Here we obtain the meson masses and the wave function which is used to smear 
the source and sink operators. 



Smearing of mesonic operators 

To fix the exponential smearing function we measure the wave function at zero temperature, 
given by ( |2l| ) with point smearing u) p (y) = 8(y) both for the quark and anti-quark at the source: 
and define 

u e {y) = G$\y,t)/G$\0,t)\ t »i = exp(-aif), (35) 

where a, p are fitting parameters. The fitted values of a and p of the Ps meson wave function 
are listed in Table g. 

To extract the effective mass from the correlators following eq. (28) we have several possi- 
bilities depending on the choice of mesonic operators $j^(t). We call m^"^ (t) the effective 
mass extracted from correlators smeared with lo and lo' at the source and the sink respec- 
tively. In figure [l] we display the effective masses for n a = 0.086. In all cases, the sink lo' is a 
point-point operator and we show three choices of lo at the source: point-point, point-exp and 
exp-exp (pp, pe and ee respectively in what follows). For S and A channels, only the correlators 
with exp-exp source are shown, since with other sources statistical fluctuations are so large that 
no clear plateaus are found. It is clear from the figure that the effective masses from different 
operators converge to the same value at large t, the worst behavior being observed for the non- 
smeared pp-pp correlator. Considerable improvement is observed for the masses extracted from 
smeared operators. Of them m^^ pp \t) is the one that converges most rapidly to a plateau, 
though it increases slightly at early stages, which is due to the fact that the source is slightly 
too wide, as discussed in section ^L^. The amount of optimization achieved by the exponentially 
smeared operator has to be analyzed through the diagonal correlator ($^(t)<I>^(0)) which is 
a sum of positive contributions from the different states - see (f^) , (||) . In Fig. g we show, in 
addition to the off-diagonal ee-pp effective mass, the masses extracted from diagonal ee-ee and 
pe-pe correlators. They both show a very similar behavior and do not reach a plateau up to 
t ~ 10, where they merge with the ee-pp result. This is an indication that the "good" behavior 
of {t) is partly due to an "accidental" cancellation of contributions from higher excited 

states which in a non-diagonal correlator may have alternating signs. Therefore also the m e ^{t) 
extracted from such correlators is no longer bounded from below by the meson mass. 

We have tried to improve the mesonic operator by performing a variational analysis in the 
basis of operators {$ pp , $ pe , $ ee }. This amounts to a diagonalization providing the best ap- 
proximation to the ground and two first excited states within the space of operators we have 
used. The result of the diagonalization is also shown on Fig. ^ where the effective masses of 
the "ground" and the "first excited" states are displayed (the "second excited" state suffers 
from large fluctuations). As can be seen, within this basis of operators no improvement is 
obtained. The ground state effective mass is very similar to the ee-ee and pe-pe ones, prob- 
ably an indication that the basis of operators used is too correlated to provide any further 
improvement. 
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IF 


a 


P 


m Ps 


my 


m s 


m A 


0.081 
0.084 
0.086 


4.05 
3.89 
3.78 


0.3785(33) 
0.3797(31) 
0.3800(25) 


1.289(8) 
1.277(8) 
1.263(8) 


0.1777( 8) 
0.1493( 9) 
0.1341(10) 


0.1962(10) 
0.1747(12) 
0.1644(13) 
0.1225(16) 


0.302(5) 
0.285(7) 
0.280(9) 
0.248(13) 


0.314(5) 
0.300(6) 
0.296(6) 
0.269(9) 



Table 2: The meson masses at zero temperature in units a~ . For the scalar channel, only 
the connected part is evaluated, a and p are the fitted parameters of the observed Ps wave 
function, and they are used to smear the quark source. 

Summarizing the observations of the source dependence and the variational analysis at 
T = 0, within the basis of three operators we have at present, the analysis does not achieve 
further optimization of the correlators. Since the effective masses extracted from all sources 
approach the ee-pp one and the latter reaches earlier a plateau we shall use the ee-pp correlator 
for the coming discussion. One should however keep in mind that such correlator is not really 
optimized in the sense of been constructed from a sufficiently optimized meson source at T = 0. 
There is a priori no guarantee that the cancellation taking place at T = will still remain at 
T > 0. We will use the departure of m^ - * from flatness as an indication that temperature 
effects start to become relevant. To control the uncertainties introduced by this choice we keep 
measuring correlators with the different types of operators, to investigate the effective masses 
source dependence also at T > 0. As long as the effective masses extracted from different 
sources converge to the same value, temperature effects will be small and the extraction of the 
meson mass from the ee-pp correlator will be safe. When the source dependence starts to be 
important we will rely in addition on other type of analysis, like the study of the wave function 
and the comparison with the free quark case. 

Spectroscopy at T=0 

We briefly summarize here the meson spectroscopy. We extract the masses at T = from the 
ee-pp propagators for simplicity (see the discussion above). For the pseudo-scalar and vector 
meson, the masses are extracted from a fit to a single exponential in the range t = 27 — 36 
(where the three types of correlators coincide; although the ee-pp propagator reaches a plateau 
much earlier, precision was not lost by this limitation). For the S and A channels, the statistical 
errors are much larger than for Ps and V, therefore we adopted the fitting region t = 12 — 20 
for the former. As was already mentioned, only the connected part of the scalar channel is 
evaluated, hence the result is only useful for a comparison with the finite temperature results. 
The values obtained are listed in Table ^. These values are consistent with the result of the 
variational analysis, where we extract masses from the fitting region t = 12 — 16. 

Masses are extrapolated to the chiral limit with the definition of k in eq. (|33"|), First, the 
pseudo-scalar meson mass squared is extrapolated linearly in 1/k to determine 1/k c at which 
the Ps meson mass vanishes. This gives k c = 0.17144(11). Then the other meson masses are 
extrapolated linearly to 1/k c . The results of the extrapolation are also listed in Table 

Comparing with the physical value of the p meson mass, the vector meson mass at k c defines 
the lattice cutoff as a~ l {rn p ) ~ 5.7 GeV. This value is about 27 % larger than a~ l = 4.5 GeV 
from the string tension. The discrepancy is consistent with results on isotropic lattices, and is 
mainly explained as an 0{a) effect. 
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Figure 1: The effective masses of correlators with various source smearing functions and the 
point sink for n a = 0.086, at N T = 72. 
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Figure 2: Effective masses of Ps and V correlators for n a = 0.086 at iV T = 72. In the left plots 
for each figure, propagators with various source and sink smearings are compared to observe 
the dependence on smearing. The right plots show the result of the variational analysis: the 
effective masses of the "ground" and the "first excited" states are displayed. 

5 Non-zero temperature 

In this section, we study how the temperature changes the meson correlators. First we shall 
compare propagators with various source and sink operators. Then we shall discuss in more 
detail what we can learn from the temperature behavior of the effective masses. For a further 
insight in the temperature effects on the meson correlators we study the t-dependence of the 
wave functions. As a result we find that the two quarks tend to stay together even in the 
deconfining phase (at least for Euclidean time scales ~ ^/T). Finally we study the temperature 
behavior of temporal ("pole") and spatial ("screening") masses which could be associated with 
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the putative quasi-particle (resonance?) states suggested by the first steps of the analysis. 

To disentangle perturbative from non-perturbative effects we shall repeatedly compare the 
measured ( "full" ) meson correlators with "free" meson correlators made out of unbound quark 
propagators. For the latter we just use a free quark ansatz and only allow the quark mass to vary 
with the temperature. This means that we consider quark-antiquark correlation functions in the 
corresponding mesonic channels in lowest order perturbation theory but with a temperature 
dependent quark-mass. As has been shown by a more involved analysis, including further 
thermal effects in the HTL approximation does not significantly change the result 24]. 



5.1 Source dependence of the propagators at non-zero temperature 
N T = 20 (below T c ) 

Let us start with N T = 20. The system is in the confining phase, hence we expect the hadronic 
spectral function to still have narrow peaks corresponding to the bound states. Indeed the 
situation is very similar to the T = case. Fig. || compares the effective masses with several 
choices of source and sink smearing functions (we have also measured here the effective mass 
with wall source). For the Ps and V channels, the ee-pp correlator appears most flat showing, 
as for T = 0, a rather clear plateau. The two diagonal correlators, smeared both at the source 
and the sink, have strongest contributions from excited states but merge with the ee-pp result 
at about t ~ 7. (Since the smeared sink suffers from large fluctuations the discrepancy of the 
effective masses at t = 9 of the sink smeared propagators and the propagators with point sink is 
probably the result of insufficient statistics.) Here again a variational analysis does not provide 
any improvement on the ee-ee and pe-pe results. 

The convergence of the effective masses from non-diagonal correlators with point-point sink 
has not taken place yet at t ~ N-t/2. This could in principle be a first signal of temperature 
effects but notice that the difference of effective masses at this value of t is of the same order 
as that observed at T = at the same t slice. 

In the S and A channels the statistical fluctuation of the effective masses are much larger 
than for Ps and V. In these channels, the effective masses from wall-pp and ee-pp correlators 
coincide. Again we observed large fluctuations at large t for the sink smeared propagators. 

Like for T = 0, at N T = 20 we conclude that the most reliable estimate of the meson 
masses is again obtained from the ee-pp propagators. The extracted masses are discussed in 



subsection 5.4 



iV T = 16 and 12 (above T c ) 

At N T = 16 and N T = 12, the propagators are very similar between the different channels 
which is a clear indication of chiral symmetry restoration immediately above T c (see section 



5.4 for further discussion of this point). Figure [4| shows the effective mass dependence on the 
mesonic operator at N T = 16 and N T = 12 for the Ps channel (the effective masses in the 
other channels are similar). To see the quark mass dependence, results for two values of K a are 
shown. An interesting feature is that the K a dependence is small. 

Although above T c the ee-pp effective masses show no longer such a clear plateau (further 
discussion on this point will be done in section 5.2) they still seem to merge with those coming 
from diagonal ee-ee and pe-pe correlators at about t ~ 4 (statistical fluctuations do, however, 
not allow a quantitative estimate). There is, however, a clear difference in the behavior of the 
non-diagonal effective masses here, as compared to T < T c . In particular the wall-pp mass 
looks as flat as the ee-pp but seems to converge to a rather different value (the difference at 
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Figure 3: Source dependence of effective masses at N T = 20 for K a = 0.086, all channels. 
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Figure 4: Source dependence of Ps effective masses at N T = 16 and 12 for n a = 0.081, 0.086. 



t ~ Nf/2 being here considerably larger that the corresponding one at the same time slice for 
N t = 20). 

The observed stronger source dependence above T c might just be an effect of periodicity or 
contamination from excited states but it is peculiar that this behavior sets in precisely at T c . 
In view of the discussion of section 2.1, we consider this as indication that strong temperature 
effects develop above T c (further discussion on this point will follow later). We have also 
performed at this Nt's a variational analysis which again just reproduced the values of the 
ee-ee and pe-pe correlators. 



5.2 Effective mass of ee-pp propagators at nonzero T 

We shall first discuss the T— dependence of propagators. Figure [| shows the effective masses 
of Ps and V meson correlators with ee-pp smearing at N T = 20, 16 and 12 for n a = 0.086. 
While at iV T = 20 a plateau appears, at N T = 16 and 12 the effective masses are no longer 
flat, their t-dependence increasing with the temperature. This holds already at t-values at 
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Figure 5: Effective masses of the correlators with exp-exp source and point sink at finite 
temperature for k g = 0.086. Shown are also the effective masses of the "free" mesons, composed 
of the free quark propagators. 



which the N T = 20 data have clearly reached a plateau (t — 4) and indicates that we have 
to do here with strong temperature effects. Due to the uncertainties related with the ee-pp 
smearing it is not possible to say whether this behavior means that the mesonic state has 
become metastable above T c , or it has been replaced by a collective excitation of increasing 
width, or that the effects of the contamination with other states from the insufficiently tuned 
source have become very strong. Nevertheless it is remarkable that we find a clear change in 
behavior setting in at T c , although smoothly connecting to the behavior below T c (at least for 
the Ps and V propagators: the S and A correlators change more significantly, in accordance 
with the chiral symmetry restoration). The same signal is provided by the observation of the 
source dependence (see previous section). Notice that the contamination with other states 
due to the imperfection of the source alone would be expected to produce a rather continuous 
dependence on the temperature, and not the clear difference in behavior observed below and 
above T c . 

Since above T c we are in the QGP phase a first thing to test is to look for signals of 
the perturbative, high T regime, that is for unbound quark - antiquark propagation in the 
mesonic channels. On figure [5] we show effective masses calculated from "free" meson correlators 
constructed out of free quark propagators where only the mass is assumed to vary with the 
temperature (as already noticed, a consistent HTL calculation [24] does not significantly affect 



the result). The free quark propagators are calculated with 7^ = 5.3 and using the same 
exponential source as for the genuine mesons. The plotted results for the "free quark mesons" 
correspond to the assumption of a thermal mass for the quarks: 



merm !? T _ . 36— ar 1 (36) 
q V6 T c 



therm 



tuned such as to give a good agreement in the p channel above T c with the measured ("full") 
effective masses. We observe that above the critical temperature the measured Ps and V effec- 
tive masses change their order and increase oc T, feature shared with the "free" mesons. The 
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inversion in the order of the masses essentially implies that the propagators are not dominated 
by a narrow state (compare |25[|), however not every wide spectral function leads to this inver- 
sion, therefore this similarity may have significance. After we tuned the free quarks to simulate 
the (full) data in the p channel, the pseudo-scalar remains, however, still well below the "free" 
results (a similar observation has been made in 0]). 

Also the stronger source dependence observed above T c is a feature which the measured 
("full") meson propagators share with the free ones - but for the latter this is much more 
pronounced and for particular sources quite different from the full mesons. For instance, for 
the wall source, the effective mass of "free" mesons is completely flat and its value is twice the 
quark mass value, while the full meson effective masses are significantly larger and vary with t. 

In the analysis of the propagators we observe therefore competing features, which hint to 
some contributions from "unbound" quarks but cannot be explained only in terms of the latter. 
To investigate further this problem we analyze the t-dependence of the wave function in the 
next subsection. 

5.3 Wave functions 

We consider the normalized wave function at the spatial origin: 

^\x,t) = G^(x,t)/G^(6,t). (37) 

We recall that these correlators are obtained in the Coulomb gauge. If the correlator (x, t) 
is dominated by a bound state <p(x, t) should stabilize with increasing t, approaching a certain 
shape. If there is no spatial correlation among the quarks, in particular in the case of a "free" 
correlator (constructed from free quarks), the corresponding cp(x,t) should become broader in 
x with t (or at best reproduce the source, for m q — ► oo) Q We shall use different sources, u, 
and observe the evolution of (p^\x,t) with t. 

Figure |6] shows the change with t of ip( ee "> (x, t) for the Ps correlator with exp-exp source at 
N T = 20 and 12, and K a = 0.086. The "free" meson "wave function" (p\^ e \x,i) is shown for 
comparison. We see now a very clear difference: The normalized wave function of the "free" 
meson shows the expected behavior, becoming increasingly broad with t at all temperatures. 
The measured wave function, on the other hand, shrinks very fast from the (wider) distribution 
implied by the exp-exp source at t = R and stabilizes very early to a well defined shape with 
increasing t. Remarkably enough, this behavior is not only seen at 0.93 T c (iV T = 20) - where, 
as expected, the wave function is similar to that at T = — , but also at 1.5 T c (N T = 12): the 
wave function above T c is only slightly wider than below. 

In Figure ^ we plot ip^(x, t) vs t for various sources, at distance x = (3,0,0) and x = 
(5,0,0). Again we show both the measured correlators and the "free" ones. As noticed, a 
genuine wave function would be represented by ratios (f(x, t) independent on t for large enough 
t, smaller than 1 and decreasing with r = \x\. 

5 For a simple illustration consider two nonrelativistic quarks of mass m a , mt and individual initial Gaussian 
distributions ip q (y,0) tx exp(— y 2 /2a 2 ) and ipq{y,Q) oc exp( — y 2 /2b 2 ), respectively, then it is easy to see that the 
width of the distribution in the relative distance x, ip(x, t) = J dy dz S(z — y — x) ip q {y, t) ipq(z, t) (essentially, our 
tfi(x,t)) develops as T(t) 2 = T(0) 2 + (l/m a + l/m b ) t , with T(0) 2 = a 2 + b 2 . Nodes in the original distribution 
may lead to occasional cancellations, but the general picture is the same. It would be physically quite unplausible 
that uncorrelated propagating quarks would tend toward a distribution in the relative coordinate narrower than 
the one they start with, whatever their spectral functions might be. 

6 The distance distribution of the quarks at t = differs somewhat from the smearing function of the corre- 
sponding source, this does not modify, however, the general features. 
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Figure 6: i-dependence of the wave function <p( ee )(x,i) at N T = 20 and 12 for n a = 0.086, 
Ps channel. The "free" wave function is also displayed. Both measured ("full") and "free" 
correlators are smeared with exp-exp functions. We represent by crosses the convolution eq. 
( |22| ) which gives the separation distribution in the source. 
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Figure 7: Ps normalized wave function ^(x,t) for x = (3,0,0) (left plot) and x = (5,0,0) 
(right plot) vs t at A^ T = 20 and 12 for K a = 0.086 and for various sources (point-point, exp-exp 
and wall). The "free" wave function is also displayed (for wall this is simply 1, independently 
of x and t). 

From the figures we see that at T = 0.93T C the measured wave function approaches with 
increasing t a unique shape, independently of the source. At T = 1.5T C the exp-exp source 
also appears to stabilize, while the point-exp and the wall sources, although indicating some 
tendency towards the same shape, do not show clear stabilization. We see therefore here a 
more pronounced source dependence, very similar to what happened with the effective masses 
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above T c . Although the tendency to an increased source dependence is a feature remembering 
of the "wave function" of the free quark mesons, at all T the measured wave function is very 
different from that given by free quarks. The latter, of course, have a broadening distribution 
in all cases, the corresponding ratios ipo(x,t) increasing steadily towards 1, and show an in- 
comparably stronger source dependence.^ The difference in behavior is particularly pregnant 
for the wall source, where the measured wave function shrinks strongly with t while the "free" 
wave function is completely flat and independent on t: ipo" all \x , t) = 1 for any m q . These 
features are observed in all four measured channels. This result is strongly indicative for the 
existence of (metastable) bound states in the mesonic channels at temperatures as high as 
1.5 T c , characterized by "wave functions" similar to those below T c . From the comparison be- 
tween the behavior at 0.93 T c and 1.5 T c we may conclude that: a) even at the highest T the 
exp-exp source still projects on a state characterized by a strong spatial correlation between 
the quarks, quite similar to the low T wave function, but b) with increasing T above T c also 
other contributions in the mesonic channels show up, without such strong spatial correlation 
(this is tentatively indicated by the increased source dependence of the wave function) . 

As already remarked, the similarity of the stabilized exp-exp wave functions seen at all T 
represents self-consistent support for our source strategy, since the latter selects a given state 
on the basis of its spatial internal structure. 

5.4 Temperature dependence of temporal and spatial masses in the mesonic 
channels 

The discussion of the previous sections has shown above T c significant temperature effects 
simultaneously with the persistence of strong binding forces between quarks. The tentative 
interpretation of these results is that even above T c (metastable) bound states are present in the 
mesonic channels. In this section, assuming the existence of such states, we try to characterize 
them by extracting from the propagators the temporal masses which would be associated with 
them (would locate the corresponding peaks in the spectral functions: "pole" masses). They are 
compared with the screening masses in the same channels and their temperature dependence 
is studied. 

Temporal masses 

As discussed in sections 4-5, at N T = 20 we extract temporal masses from the correlator 
use for computing the mass the last three points near to N T /2. The fitted values 
are listed in Table ||. Masses extracted from diagonal correlators obtained in the same fitting 
region are consistent with these values within statistical error. 

In the case of N T = 16 and 12, the situation is far less clear. Here no plateau is seen, neither 
for the diagonal nor for the ee-pp correlators. We decide to extract masses from the latter, 
again by using the last three points near to N T /2, but we should be aware that these masses 
(if they do at all characterize some states) may be misestimated. The resulting values are also 
found in Table ||. They are consistent with the result of the diagonal correlators, which suffer 
from large statistical fluctuations. We note here that we quote statistical errors only, there are 

7 We use the same free quark masses ( |36| ) which were employed for the comparison of the effective masses. 
The mass dependence is monotonous and heavy quarks just reproduce the initial distribution (e.g., the crosses 
on Fig. ^| for the exp-exp source) at all t, as expected. Of course one can find some free m q and some source 
to approach some of the data points, but not to reproduce the vast majority of the features, in particular the 
shrinking of the measured (full) wave functions with t for wide initial distributions cannot be reproduced by any 
free quark ansatz. 
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N T K a 


mp s 


my 


mp s 


my 






20 0.081 
0.084 
0.086 

K c 


0.1708(17) 
0.1455(19) 
0.1315(22) 
0.040(9) 


0.1869(17) 
0.1684(18) 
0.1595(19) 
0.1233(25) 


0.292(12) 
0.277(10) 
0.271(10) 
0.243(11) 


0.294(16) 
0.279(13) 
0.272(12) 
0.242(9) 


0.1804(18) 
0.1516(22) 
0.1354(27) 


0.2036(21) 
0.1816(26) 
0.1701(31) 
0.1267(48) 


16 0.081 
0.084 
0.086 

K c 


0.1835(14) 
0.1751(15) 
0.1723(15) 
0.1568(19) 


0.1757(14) 
0.1690(14) 
0.1678(14) 
0.1564(16) 


0.1877(13) 
0.1885(13) 
0.1888(13) 
0.1903(15) 


0.1839(26) 
0.1731(13) 
0.1747(13) 
0.1804(13) 


0.3169(34) 
0.3085(48) 
0.3036(61) 
0.2868(96) 


0.3364(27) 
0.3329(31) 
0.3312(32) 
0.3244(41) 


12 0.081 
0.084 
0.086 

K c 


0.2126(13) 
0.2100(13) 
0.2101(13) 
0.2062(14) 


0.1986(14) 
0.1979(13) 
0.1996(12) 
0.2001(12) 


0.2217(13) 
0.2255(13) 
0.2275(13) 
0.2352(13) 


0.1981(12) 
0.2032(12) 
0.2061(12) 
0.2164(13) 


0.4096(16) 
0.4036(16) 
0.3995(16) 
0.3868(16) 


0.4224(21) 
0.4180(21) 
0.4148(21) 
0.4053(21) 



Table 3: The temporal masses and the screening masses divided by £ = 5.3 at finite temper- 
ature. 



large systematic uncertainties due, among others, to the smearing function dependence of the 
correlators. Since we use a non-diagonal correlator (which does not provide an upper bound 
for the mass) we cannot say in which direction this uncertainty goes. 

In the next step the extracted masses are extrapolated to the chiral limit. At N T = 20, 
in the confined phase, we extrapolate the Ps meson mass squared linearly in k. For the other 
mesons we use linear extrapolations. The extrapolations are shown in Fig. || they are similar 
to those at N T = 72. Though the Ps meson mass at the critical hopping parameter does 
not completely vanish, this can be explained as an 0(a) error and by the uncertainty in the 
extraction of the mass. 

Above T c also the Ps meson mass is extrapolated linearly. The k dependence for all mesons 
is very small. The resulting mass values at the chiral limit are listed in Table |3|. Generally 
above T c the quark mass dependence of the meson masses is small. 



Screening masses 

Before discussing the temperature dependence of the "masses" , we briefly describe the extrac- 
tion of the spatial ("screening") masses. Since the spatial distance is large ( ~ 3 fm) we use 
for this the unsmeared pp-pp correlators and the same procedure as in the calibration (see 
the Appendix). We have verified that these masses are consistent with those obtained from 



propagators using gauge invariant smearing techniques [26|. 

We measure Ps and V meson masses at all T. At zero temperature we used a.p.b.c. in all 
directions and required the masses in the space and time directions to represent the same phys- 
ical masses. Thereafter we switched to periodic boundary conditions in the spatial directions, 
but this did not induce seizable changes in the masses. 

At N T = 20, the spatial propagators show almost the same behavior as at iV T = 72. Above 
T c , the effective masses approach a plateau earlier than below T c , the obtained values are 
therefore more reliable than in the confining phase. 

Masses are again extrapolated to the chiral limit. Also here, as for temporal masses, due to 
uncertainties and partly to 0(a) effects the naive extrapolation of the Ps meson mass squared 
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Figure 8: The chiral extrapolation of the temporal masses at N T = 20. 



does not vanish exactly at l/n c 



Temperature dependence 

The temperature dependence of temporal and screening masses is summarized in Fig. |9[ The 
values are given in units of a^ 1 . The horizontal axis is 1/N T and we have T = l/N T a T with 
a^ 1 = 4.5(2) GeV. The four points (iV T = 72, 20, 16 and 12) correspond to the temperatures 
T ~ 0, 0.93 T c , 1.15 T c and 1.5 T c . The vertical gray line roughly corresponds to the critical 
temperature. 

In the confining phase, the temporal and screening (spatial) masses coincide, however above 
T c they become increasingly different. This is to be expected, since the former (whether they 
represent bound states or not) are related to propagation in plasma with the transfer matrix 
of the original problem, while the latter correspond to a T = problem with asymmetric finite 
size effects, strongly increasing with the temperature (one of the "spatial" direction becomes 
squeezed as 1/T). In agreement with other works |12|] , screening masses above T c are cx T and 
close to twice the lowest Matsubara frequency (the lowest quark momentum in the squeezed 
direction with a.p.b.c), although remaining below it for all T. On Fig. pi] we plot the ratios m/T 
for the three non-zero temperatures together with twice the lowest Matsubara frequency for 
comparison. The temporal masses above T c are also proportional with T but with a significantly 
smaller slope. The slight decrease of the Ps and V temporal masses (xl/T) above T c in the 
upper plot of Fig. [l(] is due to the large quark mass in the simulation, which produces a term 
oc 1/T in this plot. This decrease disappears in the chiral limit (lower plot). 

As a phenomenological parameter to succinctly quantify this behavior we introduce: 
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Figure 9: Temperature dependence of masses (in a" 1 ) at n a = 0.086 (top) and in the chiral 
limit (bottom). Full (open) symbols correspond to spatial (temporal) masses. The grey vertical 
line roughly represents T c . 



Since at high T the quarks are expected to exhibit an effective (temporal) mass m e q ^ ~ g 2 T/ v6 
|27|, R ~ 1 - 0.26c/ 2 for T > T c . From Fig. [Ti we see that our data tend toward this regime 
but are still well below it even at T = 1.5 T c . Notice that because of large lattice artifacts in 
our data mass ratios are more reliable than absolute values. 

Let us note again that above T c all four channels show almost the same masses. In the 
present quenched calculation, the chiral symmetry is not involved in the dynamics and the phase 
transition is the deconfining transition. Nevertheless, the chiral symmetry seems to be restored, 
which indicates a close relation between the two viewpoints of the QCD phase transition. This 
agrees with the old observation that the chiral condensate also feels the deconfining transition 
of pure gauge theory - see, e.g., [28]. 
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Figure 10: m/T at T ~ 0.93 T c , 1.15 T c and 1.5 T c for k c = 0.086 (top) and in the chiral limit 
(bottom). Full (open) symbols correspond to spatial (temporal) masses. The grey vertical line 
roughly represents T c . The dashed horizontal line corresponds to twice the lowest Matsubara 
frequency. 



6 Conclusions and outlook 

In this quenched QCD analysis the changes of the meson properties with temperature appear 
to be small below T c . Above T c we observe apparently opposing features: On the one hand, the 
behavior of the t-propagators, in particular the change in the ordering of the mass splittings 
could be accounted for by contributions from free quark propagation in the mesonic channels, 
which would also explain the variation of rri^^{t) both with t and with the source. On the other 
hand, the behavior of the wave functions obtained from the 4-point correlators suggests that 
there can be low energy excitations in the mesonic channels above T c appearing as metastable 
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Figure 11: The phenomenological parameter R eq. (38) as function of T. The grey vertical 
line roughly represents T c . 



bound states which replace the low temperature mesons. They would be characterized by a 
mass giving the location of the corresponding bump in the spectral function. In this case the 
variation of m e jj{t) with t and with the source would indicate a resonance width for these 
states, increasing with T. Our treatment of the low energy states using the exp-exp source 
introduces, however, uncertainties which do not allow for quantitative conclusions. As we have 
seen this source is indeed slightly too wide and includes some contributions from excited states 
compensating each other in the ee-pp correlator at T < T c . At high T the "effective" mass 
becomes therefore increasingly ambiguous. Remember, however, that our source is not chosen 
arbitrarily but selects a state according to the internal structure of the latter on the basis of 
the similarity with the wave function of the T = mesonic ground state. What we find is that 
at all temperatures there is a tendency for stable spatial correlation between quarks with a 
shape similar to the T = wave function. It seems therefore reasonable to hypothesize that 
at all temperatures this source finds a low energy mode characterized by strong, stable spatial 
correlations between the quarks, and that the properties of the propagators taken with this 
source will reflect within some uncertainties the properties of this mode. 

We see clear signal for chiral symmetry restoration starting early above T c . Since this is 
a quenched calculation, this effect is completely due to the gluonic dynamics. Above T c the 
screening masses increase faster than the (temporal) masses, remaining however clearly below 
the free gas limit. The exact amount of splitting among the channels and the precise ratio 
between and may, however, be affected also by uncertainties in our £ calibration, in 
the definition of the source etc: Especially the temporal masses might be misestimated (if they 
are at all well defined). But the semi-quantitative picture of much (and, with T, increasingly) 
steeper "screening" propagators (as compared with temporal ones) is undoubtful. 

A possible physical picture is this: Mesonic excitations are present above T c (up to at least 
1.5T C ) as unstable modes (resonances), in interaction with unbound quarks and gluons. They 
may be realized as collective states, by the interaction of the original mesons with new effective 
degrees of freedom in the thermal bath, or as metastable bound states (of thermal quarks?). 
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To the extent that our results can be interpreted as supporting this picture, we should note 
that: 

- the "temporal masses" of these putative states in the pseudo-scalar and vector channels 
connect smoothly to the "pole" masses of the mesons below T c (the S and A correlators change 
more significantly, in accordance with the chiral symmetry restoration), 

- the increase of the temporal meson masses with T is to a certain extent similar to that of 
the quark thermal mass if we assume g 2 ~ 1 (while the spatial masses increase much faster), 

- the wave functions characterizing these states are very similar to those of the genuine 
mesons below T c , 

- since there is no pair creation (no dynamical quarks) a "decay channel" for these putative 
resonances may be meson — ► qth + qth where by qth + qth we indicate some strongly interacting 
but unbound quark states, 

- the incomplete (quenched) dynamics of this simulation seems to already provide enough 
interaction to account (besides for mesonic states below T c ) for chiral symmetry restoration 
and binding forces above T c ; nevertheless effects of dynamical quarks of small enough masses 
may add further features to this picture. 

Although our results are consistent with the above picture, there may be also other possi- 
bilities (cf H f|], cf and references therein). We also see agreement with the earlier study of 
meson propagators including dynamical quarks (but without wave function information) Q, 
which finds masses and (spatial) screening masses oc T above T c and indication for QGP with 
"deconfined, but strongly interacting quarks and gluons" . Investigations of mesonic correlation 
functions in HTL approximation |24j] show that the observed effects cannot be reproduced from 
perturbative thermal physics. The wave function analysis in our work indicates in fact that the 
strong interactions between the thermal gluons and quarks may even provide binding forces 
which partially correlate the latter in space. 

The complex, non-perturbative structure of QGP (already indicated by equation of state 



studies up to far above T c flCJfl , see also [29]) is thus also confirmed by our analysis of general 
mesonic correlators. From our more extended study however, especially from the, here for 
the first time investigated, spatial correlations between quarks propagating in the temporal 
direction at T > T c (wave functions), the detailed low energy structure of the mesonic channels 
appears to present further interesting, yet unsolved aspects and therefore provide an exciting 
and far from trivial picture of QGP in the region up to (at least) 1.5T C . Further work is 
needed to remove the uncertainties still affecting our analysis. This concerns particularly the 
£ calibration and the question of the definition of hadron operators at high T, which appear to 
have been the major deficiencies, besides the smaller lattices, affecting earlier results JlT[ ] . We 
are also trying to extract information directly about the spectral functions j3(| . Finally we are 
aiming at performing a full QCD analysis in the near future. 
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Appendix: Calibration of anisotropy parameters 



To determine the gauge field anisotropy £, we use the ratios of the spatial-spatial and spatial- 



temporal Wilson loops |3U 32 



Ra(r,x) = W aa {r + l,x)/W aa (r,x), 

R T (r,t) = W aT (r + l,t)/W aT (r,t). (39) 



Then the matching condition fl34|) is 

R r7 (r,x) = R T (r,t = ^x). (40) 

80 configurations are used for this analysis. 

In the determination of £, we vary the minimum value of r x x (with corresponding choice of 
t), where r < 5 (r = 1 and x = 1 are not used to avoid short distance effects). The largest value 
of x for each r is chosen with consideration to the statistical errors. We obtain £ = 5.397(22) 
(r x x > 4), 5.340(40) (r x x > 6), 5.248(101) (rxi>8) and take therefore £ = 5.3(1) in the 
following. 

To determine the lattice spacings, the heavy quark potential is measured. The extracted 
value of the string tension ^fa = 0.480(23) together with physical value ^J<J v hys = 427 MeV 
gives the cutoffs a" 1 = 0.85(3) GeV and a" 1 = 4.5(2) GeV. The spatial extent of the lattice of 
about 3 fm (~ 4 times 1/T C ) is considered sufficiently large to treat hadronic correlators. 

We then proceed to the fermionic calibration. We fix the value of n a and vary 7^ to find 
out the value which gives the same anisotropy £ as for the gauge field. We define the fermionic 
anisotropy using correlators in temporal and spatial directions, expected to behave at large 
distances like 

G { ZT PP \t) = <*£°(*) *^(0)t> -> C« exp(-mj?t) (41) 

and the same expression with time replaced by one of the spatial directions z, behaving as 
exp(—m^z) for large z. 
In the calibration, we measure the pseudo-scalar (tm = 75) and the vector (71, 72) meson 
correlators. Here we adopt anti-periodic boundary conditions in all four directions, hence the 
expected behavior is a hyperbolic cosine. The physical isotropy condition (|34]) is then applied to 
the effective masses. Figure ^ shows rn^fJt) obtained by solving (pS|) and the corresponding 
rnffJz) for k s = 0.081 and with two values of 7^. On the figure m^{z) is divided by £ = 5.3 

to be compared with mj)f(t) (i.e., it is given in units a~ l ). 

For = 4.05, the spatial effective mass divided by £ = 5.3 coincides with the temporal 
one. Although the former shows no plateau because of the small number of spatial sites, the 
temporal effective mass, which is finer spaced, does reach a plateau in the large t region. It 
is consistent to expect that if both masses agree (after £— rescaling) in the region where the 
temporal mass shows a plateau, the spatial mass is also dominated by the ground state. We 
therefore determine £p, the fermionic anisotropy, as the ratio of the spatial effective mass at 
z = 5 and the fitted value of the temporal mass in the interval t = 27-36. The value of £f 
and the extracted masses are confirmed by studying correlators with smeared operators which 
do reach plateaus much earlier. The smearing procedure of the correlators in the temporal 
direction is described in detail in the next section. For the correlators in spatial directions, 
we apply the gauge invariant smearing technique |2(| (since the configurations are fixed to the 
Coulomb gauge). 
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Figure 12: The effective masses of correlators in temporal and spatial directions for K a = 0.081 
with two values of jp, 4.02 and 4.05. 



Though these calculations are carried out with periodic boundary conditions in spatial 
directions, the dependence of masses on the kind of boundary conditions is sufficiently small 
on the present lattice. 

Figure 13 shows the dependence of £p on jp. The values of £j? from Ps and V mesons are 
slightly different, but consistent within the present accuracy. We adopt the averaged value of 
Ps and V channels and estimate the error as their difference. 

We use three sets of («;„ jp): (0.081, 4.05), (0.084, 3.89) and (0.086, 3.78). In Table § 
these values are listed together with the number of configurations used for calibration. For the 
second set, two values of jp are tried. The meson masses quoted in the table are determined 
in section 4 using smeared correlators. 

Another procedure to calibrate the fermionic action using the dispersion relation is proposed 
in |]34| . In the present case, however, the procedure used above seems more appropriate, since 
comparison of pole and screening masses at finite temperature is one of the important goals of 
this work. 
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